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Due to quantum gravity fluctuations at the Planck scale, the space-time manifold is no 

\ longer continuous, but discretized. As a result the Lorentz symmetry is broken at very 

, ^ , high energies. In this article, we study the neutrino oscillation pattern due to the Lorentz 

Invariance Violation (LIV), and compare it with the normal neutrino oscillation pattern 

due to neutrino masses. We find that at very high energies, neutrino oscillation pattern is 

^ very different from the normal one. This could provide an possibility to study the Lorentz 

O !_ Invariance Violation by measuring the oscillation pattern of very high energy neutrinos from 

D ■ 

<-j ■ a cosmological distance. 

PACS numbers: 11.30.Cp,11.30.Rd,13.15.+g,14.60.Pq 

ON 
O 

^ ; I. INTRODUCTION 

o 

Nowadays the violation of the Lorentz symmetry has attracted increasing attention, because 

the Lorentz symmetry is the main concept of special relativity and any relativistic theory which is 
• i-^ 

invariant under continuous Lorentz transformations. A growing number of speculations suggests 

H ' fill 

that Lorentz invariance might be violated or deformed at very high energies [II] PI- The l° ca l Lorentz 



symmetry has been examined in many sectors of the standard model (SM) relating to photons, 
electrons, protons, and neutrons and none of Lorentz invariance violation (LIV) has been 



identified so far in these sectors for low-energies. The Lorentz invariance should be violated at 



very high energy scale or the Planck scale, since the Lorentz group is unbounded at the high 



aoost 



(or high energy) end, in principle it might subject to modifications in the high boost limit 
The Lorentz symmetry is based on the assumption that space-time is scale-free, namely there is no 
fundamental length scale associated with the Lorentz group. However, due to violent fluctuations 
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of quantum gravity at the Planck scale M p \ = hc/X p \ ~ 10 19 GeV, A p i = ^JfiG/c 6 = 10 -33 cm, the 
space-time manifold is no longer continuous, but discretized, and as a consequence, the Lorentz 
symmetry is broken. The discretization (foam structure) of space-time manifold with a minimal 
spacing ~ A p i was first discussed by Wheeler [§], and have been intensively studied in literatures 



(see for example 



13, 0). 



In Rcf. 



Ill ] , by using the universal entropy bound, it has been shown 



that the space-time has a minimum length scale proportional to the Planck length, leading to a 
discrete space-time structure. 



The possibilities of Lorentz invariance v\o. ation have been studied in quantum-gravity models 



121 ]. string theory 



la. 1 141. Loop gravity [la. Il6j]. non-commutative geometry jl7l |- 



19[],the doubly 



special relativity (DSR) 20] . In addition, there are some other effective field theories for Lorentz 



211 ] . the minimal standard model extension 



23 



24]. 



violation, for examples, the Coleman- Glashow model 
(SME) [22J, and the newly proposed standard model supplement (SMS) 

In recent years, there has been much interest in testing LIV effects. However, observational 
tests face a major obstacle of practical nature: LIV effects due to quantum-gravity are expected to 
be extremely small because of Planck-scale suppression, and low-energy measurements are likely 
to require very high sensitivities |25(. This leads to the use of high energy astrophysics data to 



provide constraints on LIV effects. For examples, Gamma-Ray Burst (GRB) data are analyzed 



261 ] . However, high energy 



to see LIV effects on the arrival time of photons at different energies 
photons can be annihilated via pair creation with the IR background, and this limits the distances 
that high energy photons can travel, and the photon number fluxes lower for higher energies limit 

Very high energy neutrinos [28( provide an alternative to test LIV effects. Practically all current 



GRB models [29[] predict bursts of very high energy neutrinos, with energy ranging from 100 TeV 
to 10 4 TeV (and possibly up to 10 6 TeV) 30| 3l|. In addition, neutrinos with energy up to ~ 10 21 
eV are supposed to be produced by cosmological objects like GRB and Active Galactic Nuclei 
(AGN) 32|. These high-energy neutrinos from cosmological distances can open a new window on 
testing LIV effects. It was suggested that neutrinos of energies as high as 10 22 — 10 24 eV could be 
produced by topological defects like cosmic strings, necklaces and domain walls 33J. Theoretical 



34|. 



framework for Lorentz violation and neutrino oscillation probabilities is proposed in Ref. 

In this article, we formulate the discretization of space-time manifold as a hyperbolic lattice with 
the lattice spacing a, and adopt the Lorentz-symmetry-breaking Wilson operator [35(] for neutrino 



fields on the lattice. We show that neutrino oscillations depend not only on their non-vanishing 
masses, but also on the Lorentz-symmetry-breaking Wilson term in high energies. This provides 
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the possibility to test LIV effects by studying high-energy cosmic neutrinos oscillations. 

II. BOSONS AND FERMIONS ON A DISCRETE SPACE-TIME 

We first give a brief review on the energy-momentum relation of free bosons and fermions on a 
hypercubic lattice of space-time. The Klein-Gordon equation for a free boson field <f>(x, t) in 3 + 1 
dimension space-time, 

= V 2 0-m 2 ^, (1) 
where m is the boson mass. Eq. (JTJ) gives the energy-momentum relation of the free boson field 

E 2 = k 2 + to 2 . (2) 



In a hypercubic spatial lattice (time continuum), one can write [36[] 



a 2 V 2 (p -)• 0(n + a) + 0(n - a) - 2<£(n) = (d+ + d~ - 2)</>(n), (3) 

where 3-dimension vector n = a (^1,^2,^3), a is the lattice spacing and the shift operators 
d ± i^(n) = 0(n ± a). Then, the Klein-Gordon equation on the lattice is given by 

4>= (d+ + d~ -2)<f>(n) -m 2 <t), (4) 

and the energy-momentum relation on the lattice is given by 

e 2 = to 2 - 2[cos(fc ; a) " 1] , ( 5) 

Eqs. ([3]l5]) are not invariant under the Lorentz transformations. For low-energy particles k ■ a 1, 
the energy-momentum relation can approximately be written by 

E 2 = m 2 + k 2 -±-(k 4 a 2 ) + 0(k G a 4 ), (6) 

which approaches the energy-momentum relation ((2|) and the Lorentz symmetry is restored. As- 
suming a more complicate discretization of space-time, we parameterize the energy-momentum 
relation as 

E 2 = m 2 + k 2 - pk 2 (k 2 a 2 ) a + 0(k 6 a 4 ), (7) 
where the third term breaks the Lorentz symmetry. 



E* 




FIG. 1: Spectrum of the native lattice Dirac equation. 

We turn to consider the energy- momentum relation of Dirac fermions on a lattice. The Dirac 
equation in the continuum space-time 



(i <fl — m)ifj(x) = 0, 



(8) 



and the energy-momentum relation is 

E = k 2 + m 2 , -co < k < +oo, 

where k is the 3- momentum of fermions. On a spatial lattice (time continuum), one uses 

(x + a**) - tp(x - aM )] 



(9) 



36] 



2a 



(10) 



where = an^ and 



il){x + a») = ^(k)e ik » an \ 

As a result, the energy- momentum relation of fermion fields on a lattice is 

E = ± sin(fca) 
a 

For • a <C 1 the energy-momentum relation becomes 



(11) 



(12) 



3„2\ 



E = ±k + 0{k 6 a 



(13) 

the usual energy-momentum relation. However, Eq. (|12j) has a problem of fermion doubling as 
shown in Fig. 1, ka = ±7r also present fermion spices. 

A chiral symmetry breaking term is necessarily added into Hamiltonian T~i so that E vs k dose 
not have secondary minima at ka = ±7r. Considering two-component fermions on a spatial lattice 
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with Wilson term 



351 ] . the new Hamiltonian is 



H = i> ] {n)a[i){n + 1) - ip(n - 1)] 



2a 

ft 

+ 4>(n)ip(n) + — ^^(n)[2^(n) - ijj(n + 1) - ip(n - 1)], (14) 



The equation of motion for ip is 

i I? 
#(n) = -7^75[V>(™ + 1) " ^(n - 1)] + — 7o[2^(n) - i/>(n + 1) - ^(n - 1)]. (15) 

Substituting a plane wave ^ = exp(iEt — ikna) solution into Eq. (fTo]) , one obtains the energy- 
momentum relation 

^9 9 sin 2 fca , „ 9 sin 4 (A:a/2) , . 

E 2 = m 2 + ^ + 4£ 2 M-^. 16 

ar cr 

For low-energy particles /ca — >• Eq. (|16p reduces to 

^2 ^ fc 2 + m 2 + ^B 2 k A a 2 + 0{k 6 a 4 ), (17) 

where the third term (S-term) violates the Lorentz symmetry. In principle B is a free parameter 
characterizing the deviation from the Lorentz symmetry. 

Ftam gravitational theories, for low-energy particles with E « , Mpl , an energv-n—n 
relation is parametrized as [37| 

E 2 -p 2 -m 2 ~ ±£ 2 ( — — ) n , (18) 



where £2 ~ 10 9 determined by the flaring AGN 



381 ] for photons m = 0. Comparing Eq. (|17p with 



Eq. (|18p . one finds that the lattice spacing a % 10 9 /M p ^, which indicates the Lorentz symmetry 
breaking scale. We will adopt this scale to study effects of the Lorentz symmetry breaking on 
hight-energy neutrino oscillations. 



III. NEUTRINO OSCILLATIONS DUE TO LORENTZ INVARIANCE VIOLATION 

In this section we study neutrino oscillations due to the Lorentz symmetry breaking £>-term in 
Eq. (|16p . Flavor neutrinos (v e , z/„, u T ) are always produced and detected via their interacting with 
intermediate gauge bosons and in the SM. Due to the parity violation, flavor neutrinos are 

not the eigenstates of the Hamiltonian and in principle they are superpositions of the Hamiltonian 
eigenstates \vi) 

H\ui) = Ei\vi), i = 1,2,3, (19) 



where E{ are the energy eigenvalues of the type-z neutrino. Using Eq. (|17|) ultra-relativistic neu- 
trinos, the energy-momentum relation can be approximately written as 

Ei « h + ^ + ^kfa 2 + (20) 

for the type-i neutrino. 

Flavor eigenstates and Hamiltonian eigenstates (mass eigenstates) are related by an unitary 
transformation represented by a matrix U, 

3 

\u t ) = E^h), ( 21 ) 
i=i 

where the flavor index I = e,fi, r. This shows that flavor eigenstate is a mixing of the Hamiltonian 
mass eigenstates \vi), (i = 1,2,3) and vice versa. Time evolution of flavor neutrino states is given 
by 

h(t)> = e~ m %) = £ e-^UuW , ( 22 ) 

i=l 

indicating, after some time t, the evolution of these flavor neutrino states leads to flavor neutrino 
oscillations. The probability of such neutrino oscillations is given by 

P Vl ^ Vll = |(^h>| 2 = E \ U U U l'iUi*U Vj \ cos[(Ei - Ej)t + <p u ,], (23) 

where 

s (m? - m?) (B? - Bh _ _ 



and = axg^UnU^UpjUi'j) [3J|-|4l|]. In the right-handed side of Eq. ([23]) . the first term is normal 
one and the second term is due to the Lorentz symmetry breaking i?-term in Eq. (|20p . 

For a two-level system of electron and muon neutrinos (u e , z/„). The unitary matrix U is explicitly 
given by 

/ 



U 



cos 9 sin 
— sin 9 cos 



(25) 



where 9 is a mixing angle. Eq. (|2ip becomes 



|f e ) = cos + sin 9\v<i) , 

\v^) = — sm9\ui) + cos 9\v 2 ) ■ (26) 
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The Hamiltonian fj 19j) in the base of the mass eigenstates \vi) is 

m 2 /2E \ 1 / B 2 a 2 E i 
m|/2£ / 16 \ B 2 a 2 E 3 

where the leading contribution to the neutrino energy E\ is obtained by assuming pi ~ p 2 = E. 
By using Eqs. (|25|27|) the Hamiltonian in the base of flavor eigenstates is given by 





' E x | 




%mass — 






i E 2 j 





# = un mass u^ 

( 



E + m * + m2 + E\ 2 ( %±3) + + E ^l 2 



\ 



-cos 20 sin 20 
sin 26* cos 20 



AE \ 1% J \ AE 16 

where AB 2 2 = B 2 — Bf, AB 2 2 = B 2 — B\ , (B 2 > B\) and the mixing angle 9 is given by 



(28) 



tan 29 = s 2Hl \ . (29) 

U22 - Un 

IV. THE CONVERSION PROBABILITY 

Based on Eq. (123p for the system of two neutrino flavors, the conversion and the survival 
probabilities of a particular flavor of neutrino with the mixing angle 9, can be written as 

Pconv(t,U) = sin 2 29 sin 2 (|), (30) 

Psurv = 1 Pconv- (^1) 

where the oscillation phase 3> is given by 42|] 



$=[ t e(T)dT, (32) 

where tj and t are respectively the initial and final time of the evolution of the system. In the case 
for vacuum oscillations, e equals to [39|] 



In the case that the Lorentz violation is present in Eq. (|28p . e can be written as 

which shows AS 2 can also generate neutrino oscillations. The discussions and calculations are also 
applied for other two-level systems of neutrinos, (9 2 3, AB% 3 ) and (9\s, AB 2 3 ). 
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Using the scaling relation 

E = E Q {t /tf'' i = E (l + z), (35) 

where to ~ 10 18 s is the present epoch, the redshift z = (to/t) 2 ^ 3 — 1 and Eq is the energy at the 
presence epoch, z = [42]. We separate the oscillation phase (j32~1) into two parts: 

$ = &vac + *LV- (36) 

then using Eqs. (|32l34p and (|35p . we obtain the vacuum and LV phases 



3 Am 2 t 0/ | J, 
LO 

a 2 . „n o ,1 



$ mc (x,Xi) = — — — (xa -x?), (37) 



^ L v(x,Xi) = ^AB 2 E 3 t (---), (38) 
8 Xi x 

where x = t/to and X\ = ij/io- From (|37|) and (|38|) . we find that the neutrino vacuum oscillation 
does not occur, Am 2 — > 0, however neutrino oscillations due to the Lorentz symmetry violation 
take place. 

Taking Xi = 0.125, corresponding to the initial time of neutrino productions at redshift z ~ 3, 
and x = 1, we obtain 

3 Am 2 

$^(1, 0.125) ~ —— — to, (39) 
iu tiQ 

$ LV (1, 0.125) ~ -a 2 AB 2 E$t . (40) 

o 

Eqs. (|39|) and (|40|) show that for very high energy neutrinos, the LV oscillation phase becomes 
more important than the vacuum oscillation phase. 

Since neutrino detectors have a finite accuracy in the reconstruction of the neutrino e nerg y, by 



averaging Eq. (|30p over the interval AEq ~ Eq, one computes the conversion probability 42| 



I r3Eo/2 

Pconv(Eo) = — - / dE'P(E'). (41) 

AEq JEo/2 

Considering very high energy neutrinos, which are produce at z = 3, using Eqs. pi|39Ml~]) . we 
plot in Fig. (|2|) the survival probability as a function of energy Eq. We find that for large neutrino 
energies, the vacuum oscillation phase (|39|) is suppressed and its contribution to the conversion 
probability Pconvi^a vr) is almost zero, and the conversion probability P C anv{v a —> v fi) is mainly 
contributed from the neutrino oscillation phase (|40p due to the Lorentz symmetry breaking. This 
implies that any observation of high-energy neutrino oscillations indicates the Lorentz symmetry 
breaking. In addition, the neutrino oscillation pattern (E'o-dependence) due to the Lorentz symme- 
try breaking is very different from the neutrino oscillation pattern in vacuum. This might provide 
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FIG. 2: The survival probability 1 — P C onv(v a — > z/g) is plotted as a function of neutrino energy Eq, 
for different values of the Lorentz symmetry breaking scale a. Am 10- 7 eV 2 and the mixing angle 
sin 2 29 



1 



421 and IAS 2 ! w |Am 2 | 



the possibility that using high-energy cosmic neutrinos, one can study neutrino oscillation pattern 
to gain some insight into the Lorentz symmetry breaking, in connection with the study of arrival 
time delay of high-energy cosmic gamma ray due to the Lorentz symmetry breaking [37| • In addi- 
tion, from the theoretical point view, it would be interesting to see how the Lorentz violation term 
(|15H17p relate to Lorentz violation operators in effective field theories, see for example Ref. j^J]. 
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